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Abstract. The invention of laser cooling methods for neutral atoms allows 
optical and magnetic trapping of cold atomic clouds in the temperature regime 
below lmK. In the past, light-assisted cold collisions between laser cooled 
atoms have been widely studied in magneto-optical atom traps (MOTs). We 
describe here theoretical studies of dynamical interactions, specifically cold 
collisions, between atoms trapped in near-resonant, dissipative optical lattices. 
The extension of collision studies to the regime of optical lattices introduces 
several complicating factors. For the lattice studies, one has to account for 
the internal substates of atoms, position dependent matter-light coupling, and 
position dependent couplings between the atoms, in addition to the spontaneous 
decay of electronically excited atomic states. The developed one-dimensional 
quantum-mechanical model combines atomic cooling and collision dynamics in a 
single framework. The model is based on Monte Carlo wave-function simulations 
and is applied when the lattice-creating lasers have frequencies both below (red- 
detuned lattice) and above (blue-detuned lattice) the atomic resonance frequency. 
It turns out, that the radiative heating mechanism affects the dynamics of atomic 
cloud in a red-detuned lattice in a way that is not directly expected from the 
MOT studies. The optical lattice and position dependent light-matter coupling 
introduces selectivity of collision partners. The atoms, which are most mobile 
and energetic, are strongly favored to participate in collisions, and are more often 
ejected from the lattice, than the slow ones in the laser parameter region selected 
for study. Consequently, the atoms remaining in the lattice have a smaller average 
kinetic energy per atom than in the case of non-interacting atoms. For blue- 
detuned lattices, we study how optical shielding emerges as a natural part of the 
lattice and look for ways to optimize the effect. We find that the cooling and 
shielding dynamics do not mix and it is possible to achieve efficient shielding with 
a very simple arrangement. The simulations are computationally very demanding 
and would obviously benefit from the simplification schemes. We present some 
steps to this direction by showing how it is possible to calculate collision rates in 
near-resonant lattices in a fairly simple way. The method can then be used to 
combine quantum-mechanical and semiclassical models for cold collision studies 
in optical lattices. 
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1. Introduction 

It has been known for a long time that light can exert a force on material objects. 
The first experiments which showed the effect of the radiation pressure of an 
electromagnetic field on matter predicted by Maxwell were done at the beginning of 
the last century pQ . Since then, the success of applying the light force in a controlled 
way to cool gaseous atomic clouds to temperatures around, and even below, the [iK 
range has given a huge impetus to the field of cold atomic physics. 

Laser cooling and trapping of neutral atoms has been a very rapidly developing 
field of physics since the mid-eighties when experimentalists succeeded for the first 
time to optically trap laser cooled atoms More recent milestones include 

the experimental realizations of atomic Bose-Einstein condensate (BEC) [3], Fermi 
degenerate dilute quantum gas HJ, superfluid-Mott-insulator phase transition in 
optical lattices |SJ, molecular BEC 0, and evidence for superfluidity in an atomic 
Fermi gas jjj. 

In general, the slow motion of the laser cooled atoms changes drastically 
the collision dynamics of the atoms compared to room temperature gases. The 
collisions may become inelastic and the research on cold collisions in magneto-optical 
traps (MOTs) has shown how the collisions consequently limit the densities and 
temperatures of atomic gases in MOTs [HIE]. The purpose of this article is to describe 
and review the work which has been done to extend the regime of cold collision studies 
to optical lattices. 

This tutorial gives a simple introduction to one specific cooling and trapping 
scheme for neutral atoms: near-resonant, dissipative optical lattices, and concentrates 
on the description of light-assisted cold collisions in this system. A thorough 
introduction into the field of laser cooling, trapping, and optical lattices can be found 
from a text book [TO|, various review articles [HllISIEllIllESIIElinillEI, and 
summer school lecture notes ^HJOH^I- Cold collision theories and experiments are 
reviewed in detail in Refs. [El HUES- 

The article is organized as follows. Section [3 gives a short introduction to near- 
resonant, dissipative optical lattices. We describe both the red- and blue-detuned 
lattices. The latter are sometimes in the literature referred also as gray or dark 
optical lattices due to the reduced number of scattered photons compared to " bright" 
red-detuned optical lattices. Section |3| describes the basic cold collision mechanisms 
in the presence of a near-resonant light, radiative heating by the red-detuned light 
and optical shielding by the blue-detuned light. Section 01 shows a formulation of the 
cold collision problem in optical lattices and overviews the central results. Finally, we 
conclude with a few remarks in Section 

2. Optical lattices 

A disordered cold atomic gas can be arranged into an ordered structure by periodic 
optical potentials which are created with suitably chosen set of interfering laser beams 
[121 1141 IT51 1161 IT71 118) . An optical lattice formed by optical potentials does not only 
trap the atoms, but it is also able to cool them via optical pumping and Sisyphus 
mechanism |2HI I24| . 

The typical temperature range for alkali metal atoms cooled and trapped in 
optical lattices falls between the so called Doppler limit Td and recoil limit T r . The 
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Doppler limit gives lower bound for Doppler cooling, the main cooling mechanism in 
MOTs, and can be written as 

T d = hT/2k B , (1) 

where fi is Planck's constant h divided by 2ir, T the atomic linewidth, and ks 
Boltzmann's constant. The recoil limit in turn describes the lower bound for 
polarization gradient cooling mechanisms, and T r corresponds to the amount of energy 
of single photon absorption or emission recoil, and is given by 

T r = (hk r f /Mk B . (2) 

Here k r is the wavenumber of the cooling lasers, and M the mass of an atom. The 
corresponding energy unit 

E r = {Tik r f/2M, (3) 

is consequently called recoil energy. A typical alkali element used for laser cooling, 
Cs, has T d = 120^K and T r = 0.2/xK. 

The modern work on optical lattices was preceded by the proposal of Letokhov to 
trap cold atoms in one dimension by using a standing light wave |25) . and accompanied 
by the experiment of Burns et al. who created crystal structures of microscopic 
dielectric objects suspended in liquid by standing light waves j^Hl- Because of the 
purity and ease of control in optical lattices, the trapped atoms are almost ideal 
for the study of phenomena which are familiar from solid-state physics. A few 
examples of these phenomena, which have been realized in experiments, include the 
Bragg scattering [23 [2H[, Bloch oscillations [221, and Wannier-Stark ladders [371] . 
Other interesting observations include the quantum Zcno effect [3JJ and dynamical 
tunneling |32| . During the last few years, the amount of research exploiting optical 
lattices has hugely increased. This is mainly because now it is possible to load far- 
off resonant optical lattices with Bose-Einstein condensates (BECs). This has led to 
the observation of superfluid-Mott-insulator phase transition in far-detuned optical 
lattices jS] which paved the way e.g. to the creation of multiparticle entanglement by 
exploiting controlled collisions [33]. For dissipative optical lattices, a so called double 
lattice creates interesting prospects for the future studies pH] . In a double lattice, the 
atoms on two ground hyperfine states can be trapped and moved in a controllable way 
with respect to each other. This opens a new avenue to study inelastic cold collisions 
in the presence of near-resonant light in dissipative lattices. 

2.1. Red- detuned dissipative optical lattices 

The simplest case of cooling in optical lattices can be described by using the atomic 
level structure with the ground state angular momentum J g ~ 1/2 and the excited 
state J e = 3/2 |23J. A single atom has two ground state sublevels \g±\/2 > and 
four excited state sublevels \e± 3 / 2 > and \e±i/ 2 >, where the half-integer subscripts 
indicate the quantum number m of the angular momentum along the z direction. 
This is shown in Fig.^with the appropriate squares of the Clebsch-Gordan coefficients 
which describe the relative strengths of the light-induced couplings between the various 
levels. 
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Figure 1. The level structure of a single atom for a red-detuned lattice. The 
squares of the Clebsch-Gordan coefficients for various transitions are shown, and 
the laser frequency uil is detuned a few atomic linewidths below the atomic 
resonance frequency loq for efficient Sisyphus cooling. The detuning of the laser 
is described by S = u>l — <+>o ■ 



2.1.1. Sisyphus cooling The laser field consists of two counter-propagating beams 
along the z-axis with orthogonal linear polarizations and with frequency ujl- The total 
field has a polarization gradient in one dimension and reads (with suitable choices of 
phases of the beams and origin of the coordinate system) 



E(z, t) = £ (e x e ikrZ ~ te y e-^ z )e-^ Lt + c.c, (4) 

where £q is the amplitude and k r the wavenumber. With this field, the polarization 
changes from circular o~~ to linear and back to circular in the opposite direction a + 
when z changes by Al/4 where Al is the wavelength of the lasers. The periodic 
polarization gradient of the laser field is reflected in the periodic light shifts, i.e., AC- 
Stark shifts, of the atomic sublevels creating the optical lattice structure. Figure El 
displays a schematic view of the optical potentials and shows how the lattice wells 
coincide with the points of pure circular polarization of the light field. 

When the atomic motion occurs in a suitable velocity range, the optical pumping 
of the atom between the ground state sublevels reduces the kinetic energy of the 
atom 23 . After several cooling cycles the atom localizes into an optical potential 
well, i.e., into an optical lattice site. Figure |3| shows the optical pumping cycles 
between the ground state sublevels and the oscillations of the atomic wave packet 
after localization into a lattice site. 

The intensity of the laser field and the strength of the coupling between the field 
and the atom is described by the Rabi frequency 

Q = 2dS /h, (5) 

where d is the atomic dipole moment of the strongest transition between the ground 
and excited states. The detuning of the laser field from the atomic resonance is given 

by 

6 = uj l - ujq, (6) 

where ojq is the atomic resonance frequency. As a unit for Q and 5 the atomic linewidth 
r is commonly used. 

The Hamiltonian for a single atom moving in the laser field given in Eq.Q is 
after the rotating wave approximation 

„2 
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Figure 2. Schematic view of the optical potentials for the two ground state 
Zeeman sublevels in a red-detuned lattice with the atomic level structure J g = 
1/2, J e = 3/2. The lattice structure is created due to the periodic polarization 
gradient of the laser field, and the points of pure circular polarization are indicated 
by dotted lines. 



Here, p 2 /2M is the kinetic energy, S the detuning of the laser, P e — Y^m=-3/2 l e "^ ( e ™l 
is the projection operator onto the excited state, and the interaction between a single 
atom and the field is 

V = - i^sm{k r z) ||e 3 / 2 )(Si/2| + ^|ei/ 2 ) (.9-1/2I 

+ ^cos(k r z) j|e_ 3/2 )(.g-i/ 2 | + -^=|e_ 1/2 )( ffl/2 || + H.c, (8) 

where z is the position operator of the atom. The spatially modulated optical 
potentials read 

U- = U sin 2 (k r z), 

U+ = U Q cos 2 (k r z), (9) 

for ground states m g — —1/2 and m g = +1/2, respectively |33]. Here, Uo is the 
modulation depth of the lattice. 



2.1.2. Localization in lattice When the steady state is reached after a certain period 
of cooling, atoms are, to a large extent, localized into the lattice sites. The optical 
lattice is in the oscillating regime if the atom completes on average more than one 
oscillation in the site before being optically pumped to a neighboring site. For less 
than one average oscillation per site, the lattice is in the jumping regime. The laser 
parameters f2 and 8 determine in which of the regimes the lattice is in |35j . It 
must be noted that tight localization and occupation of the lowest vibrational levels 
of a periodic lattice potential increases the optical pumping time, and the time of 
localization within a single lattice site becomes longer compared to the semiclassical 
values |23j . Since we are mainly interested in the case when the two atoms undergo 
an intra-well collision, the chosen parameters in our work correspond to the jumping 
regime of the lattice. 
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Figure 3. Sisyphus cooling and the localization of an atom into the optical 
lattice. A possible time evolution for a single atom wave packet is shown for two 
ground state Zeeman levels: (a) m g = —1/2 and (b) m g = +1/2. The result 
shows the optical pumping cycles and the localization of a single atom into the 
optical lattice. The discontinuous changes between the two ground states are due 
to quantum jump events from the excited state (not shown), selected to happen 
randomly with an appropriately weighted probability. If the run is repeated, the 
jumps would appear at different times. 



In our studies, the laser field is detuned a few atomic linewidths below [361 1371 l3"5] 
or above |39| the atomic transition. The most typical detuning range we use is 
3r < S < fOr towards the oscillating regime. Whereas the most common range 
in the dissipative lattice experiments done so far is slightly higher, 5r < <5 < 20r, and 
in the oscillating regime. 

2.2. Blue-detuned dissipative optical lattices 

Two counterpropagating orthogonally polarized blue-detuned (<5 > 0) laser beams can 
efficiently cool down atoms which have the level structure J e = J g or J e = J g — 1 |I 7 | . 
The lowest position-dependent eigenstate of the system is flat in this case and is not 
directly coupled to the light field at any point of space. Thus, despite of the cooling, 
the atoms are not efficiently trapped. 

The problem can be circumvented by the use of either transverse | l()j or 
longitudinal [H] magnetic fields with respect to the laser propagation axis. We have 
used for collision studies the proposal of Grynberg and Courtois with a magnetic field 
along the laser axis 021 > see below. We note that there exists also possibilities to 
create blue-detuned gray lattices by all-optical means 021 HH • 

2.2.1. Grynberg- Courtois gray optical lattice An applied longitudinal magnetic 
field removes the degeneracy of the atomic states of blue-detuned optical molasses 
mentioned above, and produces the necessary spatial modulation for the optical 
potentials and the lattice structure |f 7) . 

In general, the atomic states are now Zeeman shifted by the magnetic field, and 
light shifted by the laser. The ratio between the two shifts can be varied by changing 
the intensity of the laser or the strength of the magnetic field. Obviously, the two 
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Figure 4. The level structure of a single atom for a blue-detuned lattice 
with the Clebsch-Gordan coefficients of corresponding transitions. This structure 
corresponds to particular hyperfine states of 87 Rb and can be used for efficient 
cooling and trapping of atoms in the Grynberg-Courtois blue-detuned lattice in 
the antiparamagnetic regime 1421 . 



extreme regimes are: a) the Zeeman shift is small compared to the light shift b) the 
light shift is small compared to the Zeeman shift. 

In case a), where the light shift dominates, the behavior of the lattice 
is paramagnetic and the increasing magnetic field strength increases the lattice 
modulation depth. When the Zeeman shift dominates in the case b), the laser field 
produces perturbations to the Zeeman shifted states. The situation resembles now the 
traditional Sisyphus cooling scheme and the lattice is in the antiparamagnetic regime. 

We have done our collision studies in this antiparamagnetic regime. See Fig. 0] 
for the used atomic level configuration J e = J g = 1 and Fig. for the corresponding 
optical lattice structure. We label the three ground state sublevels with \g±x >, \g >, 
and the three excited state sublevels with \e±i >, |eo >, where the integer subscripts 
indicate the angular momentum projection quantum number m along the z-axis. 
Because the standing laser field has only circular components, and the Clebsch-Gordan 
coefficient between m e = and m g — states is zero, the atoms are rapidly pumped 
to the A subsystem of the whole state structure; thus, in this level configuration, 
the atoms are trapped to the ground substates which have an angular momentum 
quantum number to 9 = — 1 and m g = +1. The excited state with m e = provides a 
way for cooling optical pumping cycles between the two trapping ground substates. 

After the rotating wave approximation the Hamiltonian for the atomic system 
interacting with the laser field given in Eq. (@J is 



where z is the position operator of the atom, and the Rabi frequency £1 is defined as 




Q = 2d£ Q /V2h. 



(12) 
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Figure 5. Schematic view of the optical potentials for the two trapping ground 
state Zeeman sublevels in a blue-detuned lattice. The periodic polarization 
gradient of the laser field creates the lattice structure and the points of circular 
polarizations are indicated by <r+ and a~ . The dashed lines give the Zeeman 
shifted energy levels which the light field modifies. 

The interaction with a magnetic field in Eq. I|1U|) is 

U = ^2miMlB,\i){i\, (13) 

i 

where the sum over i includes all the ground and excited states, and the Zeeman shift 
factors Q,Bi for the two trapping ground substates m = ±1 are equal. 

The level structure which we have used in our studies for a blue-detuned lattice 
can be found in 87 Rb, which has F = 1 hyperfine states for both the 5Si/2 ground 
state and the 5P 1 / 2 excited state. This is actually the element and the level scheme 
which has been used in the blue-detuned lattice experiment of Hemmerich et al. 0U] > 
even though in their case the lattice structure and the orientation of the magnetic 
field differs from what is presented here. 

The reason for the choice we have made for the used level structure and an 
antiparamagnetic regime of the Grynberg-Courtois lattice is in their simplicity. The 
cooling mechanism resembles the traditional Sisyphus cooling, making it more relevant 
to compare the results between the red and blue-detuning studies. Moreover, it is 
necessary to use only three levels of the A-subsystem instead of all the six levels 
of a single atom. Thus the number of product state basis vectors for interaction 
studies between the atoms can be reduced in the blue-detuned case. This reduces the 
computational resource requirements and speeds up the simulations when compared 
to red-detuned case (which has to use all the six substates of a single atom for 
the interaction studies). Thus for blue-detuned lattices it is easier to make a wider 
exploration of parameter space, if required. 

2.3. Basic theoretical approaches 

Because of the polarization gradients, the laser couples the multitude of Zeeman 
substates of the atom in a position dependent way and the spontaneous emission 
caused by the coupling to the vacuum plays a crucial role in the optical pumping 
process. Thus, to describe the atomic motion in optical lattices one has to solve 
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the problem of a multi-level atom coupled to a monochromatic laser field and to a 
quantized electromagnetic environment in its vacuum state. 

It is possible to treat the external laser field classically since the fields which are 
considered weak from the lattice point of view still contain a large number of photons. 
Typical laser intensities used in experiments are a few mW/cm 2 . The treatment of the 
interaction between a classical field and an atom is typically done with the rotating 
wave approximation, which neglects the terms that do not conserve the total energy. 

The general form of the task is to solve the master equation for the density matrix 
p of the atomic system 

ih^ = [H,p}+C rel [p}, (14) 

where H is the system Hamiltonian, see Eqs. Q, and C re i includes the 

spontaneous emission part due to the coupling to the environment. 

It is extremely difficult to find the exact analytical solution for this equation even 
in the case of the simplest atomic level schemes used for optical lattices. One can 
try to find approximations which allow an analytical treatment, or the combination of 
analytical and numerical calculations to Eq. 1)14(1. Another possibility is to simulate 
the optical lattice system on a computer, especially, simulations by the Monte Carlo 
wave-function (MCWF) method may provide a convenient way to obtain the solution 
of Eq. (|14|l 45J 1461 1471 14"5] . In the case of laser cooling, the combinations of analytical 
and numerical treatments can treat 2D case 05] > but the 3D problem has been solved 
so far only by the MCWF simulations 

A common feature for analytical treatments of optical lattices is the adiabatic 
elimination of the excited states, thus reducing the number of the levels in the problem 
at hand |16|. One can then try to solve the master equation directly by numerical 
integration, introduce further approximations, or exploit the translational symmetry 
properties of the system |5T] . 

The MCWF method was originally developed for the problems in quantum optics, 
where in many of the cases the direct quantum-mechanical solution of the system 
density matrix is very difficult or impossible to obtain. The method has been applied 
e.g. to the resonance fluorescence spectrum of ID optical molasses and to 3D 
laser cooling jSHj. For more examples, see Ref. |48|. By now the method has also been 
applied outside the field of quantum optics, e.g. into transport problems in condensed 
matter physics [SSj- The key idea of MCWF method is the generation of a large 
number of single wave function realizations which include stochastic quantum jumps 
of the system studied. The final result for the system density matrix and the system 
properties can then be calculated as ensemble averages of single realizations. Section 
14.21 presents the method in more detail. 

We have chosen for our collision studies in optical lattices the MCWF method 
since it has been a widely used benchmark method for various semiclassical theories for 
cold collisions in MOTs P2j . The method gives a full quantum-mechanical description 
of the atomic system (the external laser field is still described classically) and treats 
spontaneous decay in a rigorous way. A semiclassical version of the Monte Carlo (MC) 
method has also been developed for lattice studies [HE1 and has been applied e.g. to 
study anisotropic velocity distributions in 3D dissipative optical lattices |56) . This 
variant describes the external atomic degrees of freedom classically, and has limited 
applicability when the spread of the wave packet influences the dynamics of the system. 
This is the case when an atom is tightly localized into a lattice site and the spread 



PhD Tutorial 



10 




Figure 6. Branching of a wave packet. The time evolution of the total atomic 
wave packet in position space is shown. In addition to treating spontaneous decay 
in a rigorous way, the full quantum-mechanical MCWF method can also account 
for the branching of the wave packet in optical lattices. 

of the packet affects essentially the optical pumping rate. Moreover, the semiclassical 
approach can not treat, e.g., the branching of the wave packet in optical lattice shown 
in Fig. For the branching of the packets, see e.g. discussion in Ref. |57| . 

3. Cold collisions between laser cooled atoms in the presence of 
near-resonant light 

The thermal velocities of the atoms in a laser cooled gas are on the order of centimeters 
per second, roughly four orders of magnitude less than in a room temperature gas. 
Consequently, the collision dynamics is strikingly different within the two temperature 
regimes. The kinetic energies involved in the cold collision process are on the order 
or less than the energy of the atomic linewidth. For example, in the case of cesium 
the linewidth of a typical cooling transition is 2400E r and a typical collision energy is 
around 1600E r . Slow atomic motion indicates that the decay time of the excitation 
becomes small compared to the time scale of the total collision dynamics. This allows 
new phenomena in cold collision processes to affect the thermodynamics of the cold 
atomic cloud • 

A general categorization of the dynamical interactions between atoms can be 
made by considering collisions occurring between two ground state or between one 
ground and one excited state atom. Here we describe binary collisions between atoms 
of same element, i.e., a homonuclear diatomic molecule, and the emphasis is on the 
collisions between a ground and an excited state atom which is the relevant scheme 
for dissipative optical lattices. 

In the presence of near-resonant light the excitation of a quasimolecule formed 
by the colliding atoms may occur at an extremely long range, even on the order of 
few thousands of Bohr radii a . Molecular potentials at these long ranges are usually 
labeled by the Hund's case (c) notation, where the component of the total electronic 
angular momentum along the internuclear axis is a good quantum number. The 
electron clouds of the colliding atoms do not overlap at long range, and the dominant 
interaction between the atoms is the resonant dipolc-dipole interaction. We give a 
description of the resonant dipole-dipole interaction in Section [O] earlier calculations 
for alkali-metals can be found in Refs. [^El EHj • 
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Figure 7. A schematic view of radiative heating of colliding atoms. The 
quasimolecule is excited at the Condon point r c and accelerated on the upper 
level before spontaneous decay terminates the process. The consequent kinetic 
energy gain is noted by AE 



The long range properties of molecular states and the sign of the detuning of the 
laser with respect the appropriate electronic transition define the possible consequences 
of collisions. A repulsive or an attractive character of the molecular potential at long 
range arises due to the relative orientation of the dipole moments of the colliding 
atoms. In the case of a red-detuned laser, the resonance or Condon point r c , occurs for 
an attractive state. For a blue-detuned laser r c occurs for a repulsive state. Depending 
on the magnitude of the detuning and the strength of the coupling laser, off-resonant 
excitation to a non-resonant state may also play a role |39j . 

A very active research field of its own is the photoassociation of laser cooled 
atoms to molecules. Photoassociation is typically done by using a large red detuning 
of the laser so that free atom pair is excited at r c to a well defined bound molecular 
vibrational state. Especially photoassociation in an atomic Bose-Einstein condensate 
has attracted wide interest recently |5U1 |^ We will not discuss further 

photoassociation here but the reader can find presentations of the field e.g. from 
Refs. 

3.1. Radiative heating by red- detuned light 

With red-detuned light, the population of the quasimolecule formed by the two 
collision partners may get resonantly excited into an attractive excited state at the 
Condon point r c (see Fig. [7J . The relative velocity of the atoms increases due to the 
acceleration on the attractive state until spontaneous emission terminates the process. 
When the atoms have again bounced apart due to the short range repulsion in the 
ground state, the pair may lose some of the gained kinetic energy in the reverse process, 
but to lesser degree. The overall effect is the heating of the colliding pair, and the 
escape of the atoms from the trap if the total gain in kinetic energy is large enough 

M 

Figure shows a semiclassical (SC) schematic view of the process. In some 
of the parameter regimes SC descriptions, such as the Landau-Zener level crossing 
model, can be used [S3]. When the SC models fail, full quantum- mechanical methods 
are needed. For example, MCWF simulations |641 Irj5] can be used as a benchmark 
method for simpler analytical semiclassical calculations. It should be emphasized that 
it is difficult in SC models to account for population recycling, which means that once- 
decayed population may get re-excited in strong laser fields. A comparison between 




Figure 8. An example of a radiative heating study in MOT from Ref. 1651 . 
The final momentum distributions for various initially Gaussian momentum 
distributions are shown. The initial momentum refers to the (p) of the narrow 
initial momentum distributions. The effects of radiative heating, the spreading of 
the distributions, and the momentum increase is clearly visible, especially for low 
initial relative momentum. For higher initial momentum some character of the 
initial distribution is still preserved in the post-collision distribution. 



various methods and their application range is given in Ref. |54) . 

Most of the radiative heating studies done so far have used a simple two-state 
description with one ground and one excited molecular electronic state. An example 
of MCWF simulation results for MOT from Ref. is shown in Fig. |SJ The results 
demonstrate radiative heating via the spreading of the momentum distributions. 
The effect clearly becomes stronger as the initial collision momentum, i.e., cloud 
temperature, decreases. We emphasize that our lattice studies include many attractive 
and repulsive states simultaneously, see Section IP1 

In addition to radiative heating, atoms may also escape from the trap by the 
fine-structure change mechanism. If the population survives on the excited state for 
small enough relative distance between the two atoms, the point where two fine- 
structure states have a crossing may be reached. Now, if the pair comes out from 
the collision in an energetically lower fine-structure state than the one in which they 
entered the collision, the pair gains kinetic energy by the amount corresponding to the 
fine-structure splitting of the atomic electronic states at large r. This energy difference 
is usually large compared to the trap depths, and consequently in this case the atoms 
escape from the trap [HI Ej ■ 

We deal in our studies with very small detunings, a few atomic linewidths only, 
and strong laser fields. It is therefore reasonable to assume that the effect of fine- 
structure changing collisions on the cloud temperature is very small compared to the 
effect of radiative heating processes, and we neglect the fine-structure change loss 
mechanism in our lattice studies. The assumption is also supported by the fact that 
there can be an order of magnitude difference between the fine-structure state crossing 
point and r c for the small detunings we have used. 

3.2. Optical shielding by blue- detuned light 

When blue-detuned light is used, the resonant excitation at the Condon point occurs 
to a repulsive excited quasimolecular state. This makes it possible to shield the atoms 
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Figure 9. A schematic semiclassical presentation of optical shielding. The 
quasimolccule is excited resonantly to the repulsive molecular state U e at the 
Condon point r c . Then it reaches the classical turning point rtp, and is finally 
transferred back to the ground state U g when arriving at r c again. If the transfer 
back to the ground state is not complete, the atom pair may gain kinetic energy 
as it is further accelerated by the excited state potential. In this case shielding 
is incomplete and the collision is inelastic. If the population transfer between the 
states is adiabatic, shielding is complete and the collision between the atoms is 
elastic. 

from close encounters [§], see Fig. |UJ If the shielding is efficient, collisions between 
atoms may become completely elastic. The mechanism would obviously be useful for 
preventing loss of atoms in optical lattices formed with a near-resonant blue-detuned 
light. 

In an optical shielding process, the resonantly excited quasimolecule population 
reaches the classical turning point on the repulsive excited state and the atoms begin 
to move apart again. The shielding becomes complete if all the population has been 
excited, no spontaneous decay has occurred, and all the population returns resonantly 
to the ground state at the Condon point. In this case, collisions become elastic (when 
photon recoil effects are ignored), and no heating or escape occurs due to the inelastic 
processes. Moreover, the ground state is emptied at a relatively long range, and no 
population reaches short distances where unwanted processes are possible, such as 
hyperfine state changing collisions. Thus, the possibility to use optical shielding in 
an efficient way allows to increase the occupation density of the dissipative lattice, in 
addition to the benefit of reducing the rate of scattered and reabsorbed photons and 
the consequent effects. 

In the past the MCWF simulations have described the efficiency of the shielding 
process in a MOT by a shielding measure P$ , which essentially describes the flux of the 
ground state population to the short range beyond the Condon point r c |661 1671 , 
see Fig. In the case of optical lattices, a more descriptive result is the momentum 
distribution in a steady state compared between interacting and non-interacting atoms 

Wi- 
lli addition of the MCWF simulations, Fig. displays also various semiclassical 
Landau-Zener results for the shielding measure. In the basic Landau-Zener approach 
(LZ), the calculation of the excitation probability is based on the linear level crossing 
model and assuming narrow initial wave packet in momentum space |6fil I7()| . The 
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Figure 10. An example of an optical shielding study in a MOT from Ref. 66 . 
Shown is the shielding measure Pg as the function of the Rabi frequency Q. 
Pg describes the ground state population flux at very short range. The dots are 
results of the Monte Carlo simulations and the lines presents the results of various 
semiclassical Landau- Zener approaches (see text). When the laser field is strong 
enough, the ground state is effectively emptied before the atoms have a chance to 
approach close to each other. This is demonstrated by the decreasing value of Ps 
for increasing Rabi frequencies. 



Landau-Zener model with decay (LZD) takes into account the exponential decay of 
the population from the excited electronic state whereas the Landau-Zener model 
with delayed decay (LZDD) accounts also for the re-excitation of the once decayed 
population due to the strong laser field |66j . 

In this Section we have described radiative heating and optical shielding processes 
by using the molecular state description. This is how the treatment has usually been 
done in the past when studying the heating and loss of atoms in MOTs. In our 
calculations and simulations we use the two-atom product state presentation instead. 
Our aim is to study the effect of collisions in a near-resonant optical lattice. Most of 
the photon scattering still occurs when the atoms are outside the range of a binary 
interaction. It would be an unnecessary complication to describe photon scattering 
and quantum jumps in a lattice by using a molecular basis. Moreover, we give a 
simultaneous description for laser cooling and collision processes. Hence, in our case 
the molecular basis is only used for a qualitative description of radiative heating and 
optical shielding processes, but it does not appear directly in our calculations. 

4. Collisions in optical lattices 

In the past cold binary collisions between atoms have been widely studied under 
conditions which correspond to the atoms trapped in magneto-optical traps |3 El 1221 • 
Previous cold collision research has concentrated on the effects of inelastic collisions 
on the properties of an atomic cloud, but has neglected the co-existence of the cooling 
processes. 

We are dealing with near-resonant optical lattices, and we have combined in a 
single framework the cooling and cold collision dynamics. Thus our approach includes 
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simultaneous dynamical processes of cooling, trapping, and collisions in atom gas. 

For far off-resonant lattices the possibility to control the cold collisions coherently 
has been proposed [7T] and realized in an experiment [331 . This allows the creation 
of an entanglement between the atoms, a step towards quantum computing in optical 
lattices [75]. Moreover, the observation of superfluid-Mott-insulator phase transition 
helps in filling in the sites of the far-off resonant lattice by controlled particle number 
0. For dissipative optical lattices, the creation of a double optical lattice opens 
interesting prospect for cold collision studies in the presence of near-resonant light 

El- 

Radiative heating and optical shielding studies in MOTs typically use the 
molecular state description of the binary atomic system. After choosing the specific 
excited molecular state and doing the partial wave expansion, usually only the lowest 
relative angular momentum ground and excited state are accounted. The other option 
has been to consider independent pairs of partial wave states, neglecting the coupling 
between the pairs in a weak field approximation |65| . Thus the descriptions in the past 
have been two-state models, neglecting the multitude of internal states of the atoms, 
and without the position dependent coupling between the atoms and the laser field. 
For optical lattices the internal states of the atoms have to be accounted because of 
the spatial periodicity of the coupling caused by the polarization gradients of the laser 
field. 

We use the two-atom product state basis [73] . For two six- level atoms in the red- 
detuned case this means 36 basis states. If the basis is transformed into a molecular 
one, there are manifolds of attractive and repulsive molecular states buried in our 
description. We are forced to do the calculations in one dimension since the limited 
amount of computational resources. It is also worth noting that the simple atomic level 
schemes we consider, do not allow to account for the spatially dependent adiabatic 
couplings between the atomic states which play a role in the system dynamics for more 
complex level schemes |17j . 

A resonant dipole-dipole interaction between the atoms in optical lattices has been 
studied with nondynamic approaches [7j][7J3[7i>lE7J[ZE]- Usually these studies assume 
fixed positions for the atoms and concentrate on the mean-field type descriptions of the 
lattice system. These approaches neglect the dynamical nature of the cold collisions, 
and the inelastic processes of radiative heating or incomplete optical shielding. Our 
approach includes the dynamical processes of cooling and cold collisions in the same 
framework. It is important to note that once the atoms are localized into the optical 
lattice sites, they are still able to move around in the lattice. For shallow lattices this 
is because the quantum-mechanical tunneling probability between the lattice sites is 
not negligible. For deep optical lattices, which are more relevant to our case, the 
atomic motion between the wells may be induced by the recoil effects combined with 
the optical pumping process. 

For the lattice parameters we use, we have noticed that the inter-well effects are 
negligible and our interest lies in the case when two atoms end up in the same lattice 
site and collide. The intra-well collision partners may then gain kinetic energy due 
to an inelastic collision and escape from the lattice. In the case of optical shielding, 
the possibility of making the collisions elastic and preventing the atoms from close 
encounters would also prevent the atoms to escape from the lattice. 

Numerical simulations are extremely heavy, especially in the case of red-detuned 
lattices where the level scheme can not be simplified. For the required computer 
resources, see Appendix A. We are forced to make some simplifications to our model, 
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for details see [37]. Most importantly, we have to neglect the reabsorption of the 
scattered photons. With increasing atomic density, reabsorption may heat the atomic 
cloud and cause radiation pressure to outward direction from the trap center |1(J|. 
limiting the achievable atomic densities. In the blue-detuned lattice, the number of 
scattered photons is largely reduced because the center of each lattice site corresponds 
to a completely dark point in space. For red-detuned lattices we merely describe 
the effects of collisions on the thermodynamical properties of the cloud. The full 
thermodynamics is not described since we neglect reabsorption. This poses some 
limitations for the applicability of our results in red-detuned lattices, but for the blue- 
detuned case our description is close to the complete thermodynamical description 
because the scattering is generally low. 

In the following Subsection we describe how we calculate the interaction 
matrix elements between two multistate atoms, especially the resonant dipole-dipolc 
interaction matrix elements. The whole problem is then formulated by using the 
MCWF method. For details, see Ref. [37j. The Section ends with the presentation of 
the central results. 

4-1- Resonant dipole-dipole interaction 

One of the early treatments for resonant dipole-dipole forces between two atoms was 
given already at the end of the 30's [TOj, and the retardation effects were discussed 
almost a decade later jHU]. Later on, Lenz and Meystre considered the resonant dipole- 
dipole interaction (DDI) for two two-level atoms in a standing- wave field [HI]. Our 
derivation follows their approach and makes the generalization to the multilevel atom 
case. 

The DDI is the first interaction to come into play when the colliding atoms 
approach each other in dissipative optical lattice. Since the interaction between the 
atoms is mediated by the quantized environment, the natural starting point is the 
two-atom system master equation and its damping part describing the coupling of the 
system to the electromagnetic environment [HI] [see also Eq. (|14Jl ] 

P= ~T2 I drTr f {H sf (t)H sf (T)p sf (T)-H sf (t)p sf (T)H sf (r) 
n Jo 

-H sf (r)p sf (r)H sf (t) + p sf (r)H sf (r)H sf (t)} , (15) 

where p is the reduced density matrix of the two-atom system, p s f the density 
matrix of the two-atom system and the field, H s f denotes the system-field interaction 
Hamiltonian, and Trj the trace over the field. 

We expand the electromagnetic field in the standard way 

E(r a ) =E+(r Q )+E-(r Q ) 

E+(r Q ) = ^i£(k)a k e ikr <* 

k 

E-(r a )= (E+(r a )) t , (16) 

where a k is the annihilation operator for mode k, r Q denotes the position of atom a, 
and 

£(k) = ^M Cj (k), (17) 
where ej is the polarization vector and V the quantization volume. 
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We use the center of mass and relative coordinates of the atom pair 

T3 r i+ r 2 , 1tA 
R= — - — , r = r 2 -ri, (18) 

and the notation 

S+ t q = CG q m (|e m + g )i i(g m \ + \e m + q ) 2 2(g m \) 

i n 

AS+, q = J2 CG9 md )i i(g m \ - \e m+q )2 2(g m \) , (19) 

m 

where CG^ are the appropriate Clebsch-Gordan coefficients, g is the polarization label 
in the spherical basis, and sub-indices label the two atoms. The interaction between 
the two-atom system and the vacuum electromagnetic field can now be written as 




where tj. q is the projection ej >q — €j ■ e q on the spherical basis eo, e±, and d the dipole 
moment of the atomic transition. 

One can identify the DDI interaction terms between the atoms as those having 
(n w + 1) = 1 where n w is the number of photons in the mode of the environment with 
mode frequency lu. In another words, the average photon number in the interaction 
process is zero since the DDI interaction can be viewed as an exchange of excitation 
between the two atoms via the environment vacuum field. After lengthy analytical 
calculations following [5T] , and using the arguments from Ref. |S2] , one can write down 
the expression for the three-dimensional resonant dipole-dipole interaction as 



Vdi P = - -hr < - 



2 

1 
3 
1 

72 



1 cos qor 



3 q r 
sin g T cos qo r 



1 - 2P 2 (cos0 r )] (S++S-+ + S+S— - 2S +0 S-o) 



(q r) 2 [qorf 



P 2 (cos0 r ) {S++S-++S+S— - 2S + oS_ ) 



+ 



cos g r / sin g r cos g r 
Qor V (gor) 2 (qor) 3 

P 2 1 (cos0 r ) cos0 r (S++S-0 + S+oS — - S+oS-+ + S+S-o) 



P|(cos0 r ) cos 20,. (S ++ S + S y) | 



(21) 



where go = ojq/c, P 2 is Legendre polynomial, are the associated Legendre functions, 
and r is the lincwidth of the atomic excited state. The angles r and (j> r are the 
angles of the relative coordinate r in the spherical basis. We have also introduced the 
operators 

(22) 



S+ q S^ q ' — (S +tq S_ q t + S+qS^q,) 
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Figure 11. The energy shifted ground state and the attractive excited state 
[labeled by the Hund's case (c) notation] molecular potentials of Cs2 for S = 
—3. Or. The repulsive potential manifold is not shown. 



where S™ q = 



Si 



and 



m=J g 

£ 

1=—J„ 



CG^ n \e m+q ) a a (g ri 



(23) 



Here a labels one of the two atoms. 

If the two atoms are positioned on the z-axis, the DDI potential reduces to the 
one-dimensional potential 



praxis 
dip 



-hT 



1 cos qor 
3 



sin go r cos qor 



Qor l (q r) 2 {qor) 3 

(S++S-+ + - 25 +0 5_o) ■ (24) 

It is worth noting that the interaction potential (|24|) includes the retardation effects. 
By diagonalizing Vdi P , it is possible to obtain the molecular potentials shown in Fig.lTTl 
We use the molecular basis occasionally for the qualitative description of the collision 
processes and emphasize that the calculations are done in the two-atom product state 
basis. 



4-2. Monte Carlo wave-function formulation 

We use a variant of the Monte Carlo (MC) method which was developed by Dalibard, 
Castin, and M0lmer 07J . The core idea of the Monte Carlo wave- function 

(MCWF) method is the generation of a large number of single wave function 
realizations including stochastic quantum jumps of the system studied. Quantum 
jumps occur to the available decay channels of the system whose environment is 
continuously monitored. In our case, detection of a photon corresponds to a quantum 
jump from an internal excited electronic state to the ground electronic state of an 
atom in an optical lattice. Solutions for the steady state density matrix and system 
properties can be calculated as ensemble averages of single wave- function realizations. 
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In general, we look for the solution of the master equation l|14|) whose relaxation 
part is in the so called Lindblad form [53] 

£rel [p] =-\Y, ( C i Ci P + P°i Ci ) + E °iP C l ( 25 ) 



where the summation over i include the system operators Ci and their adjoints C[ 
(we define these operators in a moment). 

To unravel the appropriate master equation by generating the realizations for the 
MC ensemble, one solves the time dependent Schrodinger equation 

iK^=H\f). (26) 
In our case \ip) describes the time-dependent two atom wave- function in position space 
Mz*,t))= J2 ^ 2 1 ;™ 1 (^,t)|jim 1 ) 1 b 2 m 2 ) 2 , (27) 

jl,j2,m 1 ,m 2 

where ji and j 2 denote the ground and excited states of atom 1 and 2 respectively, m\ 
and to 2 the z-component of the angular momentum. Due to the limited availability 
of computer resources we have to fix the position of atom 1 by setting z\ = 0, and z 2 
is then the both position of the moving atom 2 in the lattice, as well as the relative 
coordinate [3*7| . 

The non-Hermitian Hamiltonian H in Eq. i|26|) is 

H = Hs + Hdec, (28) 
where the system Hamiltonian Hs in our case includes the atom-laser interaction 
Hamiltonians expanded in the two-atom Hilbert space, and the resonant dipole-dipole 
interaction between the atoms, Eqs. I|7I1UI24|) . 

The non-Hermitian part includes the sum over the various allowed decay channels 

3i 

H DEC = -fE C l C i< ( 29 ) 

3 

where Cj are the jump operators corresponding to particular decay channels and can 
be found out from the relaxation part of the master equation l|25|l . 

During a discrete time evolution step of length St the norm of the wave function 
may shrink due to Hd ec ■ The amount of shrinking gives the probability of a quantum 
jump to occur during the short interval 5t. Based on a random number one then 
decides whether a quantum jump occurred or not. Before the next time step is taken, 
the wave function of the system is renormalized. If and when a jump occurs, one 
performs a rearrangement of the wave function components according to the jump 
operator Cj, corresponding to the decay channel j, before renormalization of \ip). 

For example, if we denote the jump of atom 1 from |e_ 1 / 2 ) 1 to | <7_ 1/2)1 as channel 
2 in our red-detuned lattice studies, the jump operator in the product state basis for 
this jump is 

C 2 = \/2/3Vf {|3_i/ 2 )i |ff-i/2)2 i(e_i/ 2 | 2(3-1/2! 

13+1/2)2 i{e_i/ 2 | 2(3+1/2! 
|e_ 3 /2)2 i(e_i/ 2 | 2(e_ 3 / 2 | 
|e-i/2>2 i(e_i/ 2 | 2<e_i/ 2 | 
|e+i/a>2 i(e-i/ 2 | 2<e+i/ 2 | 
-1/2)1 |e +3 /2)2 l(e_i/ 2 | 2(e +3 / 2 |}. (30) 
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Here, the factor y2/3 is the Clebsch-Gordan coefficient of the corresponding 
transition. After applying the jump operator Cj, the wave function is still in a 
superposition state, but it has collapsed into a subspace of the product state basis 
vectors, leaving only one ground state level component of the atom populated. 

In general, the jump probability into the decay channel j for each of the time- 
evolution step St is 

p j = st{i>\atc J \^). (31) 

Thus, the jump probability for an example channel 2 in Eq. I|30[l for each time step is 

p 2 = ^rji^T + l«l 2 + Hl--l%\ 2 

+ \^--i%\ 2 + ivc;:! 2 + ivc; 2 T} . (32) 

Reference [37J presents in detail the implementation of the MCWF method in our 
lattice studies. There is a large number of numerical problems one has to solve for 
the MC implementation of two atoms in a lattice, for a list and solutions see Ref. |37| . 
For general numerical tools, e.g. split Fourier or Crank-Nicholson methods, to solve 
time dependent Schrodinger equation, see Ref. 

4-3. Red-detuned lattices 

In this Subsection we present the main results from Refs. |36l |3"T| which deal with 
collision dynamics in red-detuned lattices. The main collision process in this case is 
radiative heating, see Section 1331 

Once the atoms are localized into the lattice sites, they are still able to move 
around in the lattice. When the occupation density of the lattice increases one can 
ask what is the effect of collisions for the cooling dynamics in optical lattices, and how 
the cold collisions affect the atomic cloud once the atoms are localized. 

At the beginning of the efficient cooling period a large fraction of atoms have 
higher kinetic energy than the optical lattice modulation depth. Atoms then have 
a high mobility and change their internal state frequently via the optical pumping 
cycles which cool them. A priori one might assume that the possible consequence of 
collisions would be a slowing of the cooling process, heating, and escape of the atoms 
from the lattice. Radiative heating studies at low temperatures in MOTs show a 
smooth widening of the momentum probability distribution corresponding to heating 
for large range of parameters see Fig. El A similar effect might be expected 

to occur in an optical lattice as well. 

It turns out that the internal structure of the atoms and the spatial dependence 
of the atom-field coupling changes the consequences of the collisions to some 
extent. Lattice structure introduces selectivity into the collision processes and atomic 
dynamics. In a lattice, the mobility of an atom between the lattice sites depends 
essentially on the kinetic energy, especially once the atoms are localized into the 
lattice sites. An atom, which has a large oscillation amplitude (corresponding to a 
large kinetic energy) in the lattice site, has a higher probability to change its internal 
ground state by optical pumping than an atom which is tightly localized into the 
vicinity of the center of the potential well. Rich dynamical features of the wave 
packet, like breathing and oscillations in a single lattice site, arise due to the fact that 
the packet here describes a superposition of the populations in the vibrational states 
of a lattice potential |84|. 
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Figure 12. An example of a single wave packet realization from the full MC 
ensemble in position space (note the direction of z-axis because of the viewing 
angle). At the beginning of the time evolution the atom oscillates in the lattice 
well which has the center at z = 0.25A. It then approaches point z = (the 
position of the fixed atom), collides with another atom, and gains enough kinetic 
energy to be ejected from the lattice. 




P(P r ) 



Figure 13. The same MC realization as in Fig. ll2l but shown here in momentum 
space. The population is transferred to the high values of momentum due to the 
collision, and the atom consequently escapes from the lattice. 

Since the high kinetic energy atoms are more mobile compared to their low kinetic 
energy partners, the high kinetic energy atoms have also a higher collision probability. 
During the same period of time the high-energy atoms change their internal ground 
state and the corresponding optical potential more often than tightly localized atoms. 
Thus, the coverage of various lattice sites, and the corresponding collision probability, 
is higher for more mobile atoms. 

With these assumptions the consequence of collisions might be simple heating of 
the atomic cloud, or escape of the atoms from the lattice. The essential ingredient 
for a large kinetic energy increase is a high excitation probability of a quasimolcculc. 
This depends on the curvature of the molecular potentials at a resonant Condon point 
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Figure 14. An example of momentum probability distributions for interacting 
and non-interacting atoms in a red-detuned optical lattice 1371 . The momentum 
p is expressed in the recoil unit p r = hk r . For interacting atoms the distribution 
gets narrower compared to non-interacting atoms. The results correspond to a 
sodium lattice with lattice the depth Uo = 339-B r , detuning <5 = — 3.01", and Rabi 
frequency Q = 2.8F 



and especially on the relative velocity between the colliding atoms. The optical lattice 
modulation depth defines the initial velocity distribution of the atoms when they 
start to move between the lattice sites after localization. It turns out that in a lattice, 
detuned a few atomic linewidths below the atomic resonance, and for lattice depths of 
a few hundred recoil energies, the surroundings are very favorable for strong excitation 
of the quasimolecule and the corresponding large kinetic energy changing collisions. 
The consequence is that the atoms mainly leave the lattice when colliding, and the 
total effect is the ejection of the hot atoms from the lattice. The ones which remain in 
a lattice have lower average kinetic energy per atom than in the low occupation density 
case of the lattice when there is no need to account for the interactions between the 
atoms. Figures El an d ED show in position and momentum spaces, respectively, an 
example of a collision which ejects the atoms from the lattice. 

The selective escape mechanism resembles evaporative cooling used to produce 
BEC in magnetic traps. Here the rethermalization of the remaining atoms is more 
limited, though. Anyhow, spatial dependence of the laser field introduces selective 
heating of the hot atoms, and the consequent escape from the lattice. It would 
probably be too far reaching to claim that this effect should be visible in an 
experiment. Our model neglects the reabsorption of photons which may affect the 
total thermodynamics of the atomic cloud at high densities, and we neglect also 
the Doppler cooling. Thus we have revealed one aspect of the thermodynamics of 
a densely-populated near-resonant optical lattice but the solution for the complete 
problem is simply out of reach for the modern computational resources. 

Figure IT~i1 shows an example of the results for collisions in a red-detuned lattice 
from Ref. [37]. This example is for a sodium lattice of depth Uq — 339 E r . The 
comparison is done between the momentum probability distributions of the interacting 
and non-interacting cases with the occupation density of 25% of the lattice. The 
central peak is clearly narrower when the interactions between the atoms have been 
included. This central peak corresponds to the atoms which are trapped in a lattice, 
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Figure 15. The momentum probability distributions for the <5 = 5r blue lattice 
case. The Rabi frequencies are: a) Q = 1.5V, b) Q = 2.0r, c) Q = 3.0r, and 
d) Q = 5. Or. The momentum p is expressed in the recoil unit p r = hk r . The 
dashed line is for the interacting and the solid line for the non— interacting atoms. 
When Q increases the regime changes from incomplete shielding, a) and b), to 
complete shielding, c), and finally to off-resonant heating in d). The momentum 
distributions also get wider due to the deepening of the lattice with increasing Q. 

and the wide wings correspond to background atoms which are presumably out of the 
recapture range of the lattice and ejected from the lattice. 

One can calculate by semi-classical means the excitation and survival probability 
for various molecular potentials. We have done a simple semi-classical analysis by using 
Landau-Zener approach in Ref. This analysis supports the conclusions presented 
above and shows the high probability for the atom pair to gain kinetic energy by the 
amount with which the collision partners are kicked out from the lattice. 

4-4- Blue- detuned lattices 

The prospect of using the trapping and cooling lasers for efficient optical shielding has 
been studied in Ref. Complete optical shielding would make collisions between 

atoms, when they end up in a same lattice site, elastic, and it would also prevent 
atoms from close encounters, reducing, e.g., inelastic hyperfine changing collisions. 
Thus efficient optical shielding could be beneficial in optical lattices in addition to 
the typical darkness of the blue-detuned lattices. The number of scattered photons in 
gray-lattices can be roughly two orders of magnitude smaller than in MOTs 00] • The 
role of the radiation pressure due to the reabsorption of photons diminishes, and our 
simplified model describes in a more realistic way the total thermodynamics of the 
atomic cloud, not only the collision aspect of the thermodynamics. 

Figure ED presents the results of the simulations where the detuning is fixed to 
5 = 5r jSHj. The results show clearly how the efficiency of the shielding changes. 
When the coupling laser is weak, a large number of the collisions are still inelastic 
ones. Wide wings appear in the momentum probability distribution, see Figs. Hoi fa) 
and II 51 (b). All the population, which is excited at the Condon point, does not return 
to the ground state when the atoms move apart again. The quasimolecule slides down 
on the tail of the repulsive state producing a mild heating effect j^E] ■ For a moderate 
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Figure 16. Momentum distributions for fixed Uq ~ 710E r . a) <5 = 1.5r b) 
6 = 5. Or c) 5 = 7. Or d) (5 = 10. Or. The off-resonant processes play a role at the 
small detuning, a). In the intermediate detuning shielding may become complete 
and the collisions between atoms are elastic, b). For larger detuning the shielding 
is incomplete, c) and d). 

field strength it is hard to see differences between the distributions for the interacting 
and non- interacting atoms, Fig. 1151 (c). The optical shielding has become complete 
and atoms collide elastically when they end up in the same lattice site. The collision 
partners begin to move apart at large internuclear distances and also the short range 
unwanted effects are avoided. A further increase in the laser field strength makes it 
possible for the population to be excited into the attractive molecular state by off- 
resonant means |67| . This is the reason for the deviation of the two distributions in 
Fig. 1151 (d), where the appearance of the wings is also qualitatively different than in 
the weak-field case of incomplete shielding. 

Figure 1161 presents another view to the shielding studies . Instead of keeping 
the laser detuning fixed, here the lattice modulation depth is kept nearly constant. It 
turns out that for a very small detuning, Fig. QUI (a), the off- resonant effects heat the 
atomic cloud. This corresponds to the regime where VL/8 > 1, and the steady state 
formation between the ground and the excited states surpasses the dynamical resonant 
excitation process. When the detuning is increased, the point of complete shielding 
is reached, Fig. El (b). In the region where fl/S -C 1, Fig. El (c) and (d), shielding 
becomes incomplete again due to the weak excitation and stimulated re-excitation in 
Condon point. 

The results demonstrate clearly that the co-existence of cooling, trapping, and 
shielding processes is possible in blue-detuned near-resonant optical lattices. The 
shielding is not always complete but by the careful choice of parameters shielding 
becomes very efficient. Moreover, this can be obtained within a typical and convenient 
parameter regime for near-resonant lattices, e.g., in Fig. El ( c ) § = 5I\ £1 = 3F, and 
Uq = 712E r , A clear advantage here is the achievement of complete shielding with 
the same lasers which provide the cooling and trapping. This is in contrast to MOTs 
where one needs to introduce additional lasers for shielding. 

Even though the available occupation densities in near-resonant optical lattices 
have been very low so far, the metastable rare-gas atoms could provide a convenient 
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case for an experimental study of shielding in optical lattices due to the clear ion signal 
that marks collision events 023 ■ 

The experimental work on optical shielding in MOTs show the saturation of 
the shielding phenomena when the intensity of the laser field is increased [§]. The 
saturation has not been present in earlier theoretical studies of shielding |661 167| . 
and we do not see the saturation of shielding here either. The results presented in 
Ref. [HHl seem to confirm the view that the saturation does not arise due to spontaneous 
emission effects [HE]. The reason for the saturation of shielding is still unclear. It has 
been attributed to various processes, in addition to the above-mentioned premature 
termination of shielding via spontaneous emission |66|. Other possibilities include 
counterintuitive or off- resonant processes involving different partial waves, or other 
processes that similarly involve multiple states (in contrast to the basic two-state 
approaches EH1 IEH EE] ) • In the case when several partial waves are considered, a 
net of Condon points (instead of only one resonant point of two-level models) opens 
a door for complicated excitation and de-excitation patterns which may affect the 
system dynamics |89|. The partial wave description is a useful frame of reference 
in normal collision studies to allow for more than one dimension. Due to the low 
collision velocities one can limit the studies to a few waves only. But the presence 
of an optical lattice forces us to abandon this description. Since our model is, due 
to the computational limitations, only one-dimensional, we can not conclusively that 
saturation of shielding should be absent in a lattice experiment. 

So far there has been very few cold collision experiments in optical lattices |9()U91j . 
These experiments showed how the lattice structure affects the transport of atoms 
by using collisions as a probe. We hope that our work serves as a motivation for 
experimentalists to do shielding studies in blue-detuned optical lattices. 

4-5. Collision rates 

The results of Refs. EH] show that in our selected parameter regime, i.e., the 
parameter regime for near-resonant lattices, the motion of the atom between the lattice 
sites (or in the lattice site occupied by the single atom only) is not strongly affected by 
other atoms, not at least for the occupation densities which we have used (maximum 
25%). Hence, binary interactions between the atoms come into play only when two 
atoms simultaneously occupy the same lattice site. This makes it possible to develop 
a method to calculate the average rate at which two atoms end up in the same site, 
and the consequent cold collision rate in a lattice, by following the trajectories of 
single atoms |2H1. Cold collision rate refers here to the average rate of atoms to reach 
the region of the resonant Condon point in the presence of near-resonant light and 
describes the rate of occurrence of radiative heating events in a red-detuned lattice, 
or optical shielding events, if blue-detuned light is used. It is assumed that the two 
atoms always collide when they end up in the same lattice site. For the parameters 
used here, this assumption is confirmed by the results in Refs. |361 137| . 

The basic idea of the developed method is as follows. A possible trajectory of an 
atom in optical lattice is given by a single MCWF realization. Trajectories in position 
space can then give information about the rate at which atoms travel over the average 
distance z a between the atoms, which in turn gives information about the binary 
collision rate in a lattice. The average distance between atoms z a corresponds to the 
mean free path of atoms between collision events in our one-dimensional model. Thus, 
if we monitor the transport of atoms in the lattice over the average distance between 
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Figure 17. The binary collision rate R for three different lattice depths Uo 
as a function of occupation density po of the lattice I38| , The points show 
the simulation results, and the solid lines the quadratic collision rate curves 
averaged from the simulation results for the specific red-detuned lattice. The 
detuning is kept constant at 8 = — 3.0T whereas the Rabi frequency increases as 
f2 = l.Or, 1.5r, 1.9r with the increasing lattice depth indicated in the legend. 

the atoms, or the atomic flux over the average distance, we also obtain information 
about the collision rates. 

Figure El displays the calculated collision rates as a function of the occupation 
density of the lattice for three different lattice depths In our simple ID case 

the collision rate does not depend on the scattering cross section and collisions are a 
measure of transport in a lattice. This was also the case in the experimental study of 
Ref. In both cases the collision rate has a quadratic behavior. 

The points in Fig. 1171 are the simulation results and the solid lines quadratic 
fits. It is interesting to note from the Monte Carlo point of view that we get results 
for a wide density range by doing simulations for only a very few values of density. 
The possibility of obtaining the result for all the values of the variable, in this case 
occupation density, from single Monte Carlo ensemble is a new feature in the MCWF 
simulations to our knowledge, at least when the MCWF method is applied to cold 
collision problems. 

Two-atom collision simulations in a lattice described in previous Subsections are 
computationally very heavy. It would be useful to find more simple means to do 
collision studies in optical lattices. The method presented above presents a step in 
this direction. For example, if the semiclassical analysis shows that for particular 
parameter values the colliding atoms have a high probability to be ejected from the 
lattice due to radiative heating, then the collision rate described here gives directly 
the loss rate of atoms from the lattice. MCWF simulations for one atom, like the 
ones reported in Ref. [2E], are fairly simple and fast to perform. This is especially true 
when compared to the two-atom case. Thus the combination of these simple one-atom 
simulations with semiclassical models for intra-well collision effects have a potential 
to simplify the studies of binary collisions in optical lattices. 
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5. Conclusions 

We have studied the cold collision dynamics between atoms in near-resonant, 
dissipative red- and blue-detuned optical lattices. The applied methods have been 
mainly based on the Monte Carlo wave-function method. A semiclassical analysis has 
been done which supports the conclusions drawn from the full quantum-mechanical 
calculations [3*?| . 

The implementation of the MCWF method to study cold collisions in optical 
lattices is not straightforward and the simulations have been very demanding from 
the computer resource point of view. This is due to the internal structure of atoms, 
coupling to the electromagnetic environment, position dependent coupling of the atoms 
to the laser field, and position dependent coupling between the atoms. 

The results for near-resonant red-detuned lattices are in quite a sharp contrast 
to the interaction studies in magneto-optical traps. Instead of a heating, a cooling 
due to the selection of collision partners from high kinetic energy atoms is seen in the 
simulation results. The blue-lattice results show the applicability of optical shielding. 
Future collision studies require simplifications, for which we propose a simple way to 
calculate the collision rate in optical lattices. 

In the past there has been many studies of cold collisions in magneto-optical 
traps, see the review jH]- The work presented in here extends the regime of cold 
collision studies into the realm of optical lattices. This is far from being a trivial 
step. The major reason is that for sub-Doppler cooling mechanisms, which exploit 
various polarization states of a laser field, it is necessary to account for the internal 
structure of atoms, and this greatly complicates the total system under study and 
the calculations. Moreover, it is not enough to formulate the problem using only the 
relative motion between the atoms in a constant laser field. The position of the atoms 
with respect to a lattice structure has to be accounted also. 

In conclusion, we have shown how research on cold collisions in the present of 
near-resonant light can be extended from magneto-optical traps to cover also optical 
lattices. In future, the recently developed double optical lattice [HI] opens interesting 
prospects for the experimental study on cold collisions in dissipative optical lattices. 
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Appendix A. Required computational resources 

The numerical simulations are demanding since we are dealing with a 36 level quantum 
system including various position dependent couplings and a dissipative coupling to 
the environment. We have used 32 processors of an SGI Origin 2000 machine, which 
has 128 MIPS R12000 processors of 1 GB memory per processor (52]. The total 
memory taken by a single simulation (fixed 5, fi, occupation density p , and atomic 
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species) is 14 GB, and generating a single history requires 6 hours ol CPU time in 
red-detuned lattice studies. A simulation ol 128 ensemble members then requires a 
total CPU time which is roughly equal to one month. The normal clock time is, ol 
course, much shorter (roughly 22 hours) since we take advantage ol powerful parallel 
processing lor which the MCWF simulations suit very well. 
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